Polymer dynamics, fluorescence correlation spectroscopy, and the limits of optical 
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In recent years, fluorescence correlation spectroscopy has been increasingly applied for the study 
of polymer dynamics on the nanometer scale. The core idea is to extract, from a measured auto- 
correlation curve, an effective mean-square displacement function that contains information about 
the underlying conformational dynamics. The paper presents a fundamental study of the applica- 
bility of fluorescence correlation spectroscopy for the investigation of nanoscale conformational and 
diffusional dynamics. We find that fluorescence correlation spectroscopy cannot reliably elucidate 
processes on length scales much smaller than the resolution limit of the optics used and that its 
improper use can yield spurious results for the observed dynamics. 

PACS numbers: 87.64.-t, 87.64.kv, 87.57.cf 

Keywords: fluorescence correlation spectroscopy, polymer dynamics, Zimm dynamics. Rouse dynamics 



en 



Oh 



0) 
O 

O 

Oh 



> 

o 
en 

CO 

o 

(N 



% 



The advent of single molecule fluorescence techniques 
has allowed for the direct visualization and measurement 
of the conformation and dynamics of individual polymer 
molecules, see e.g. [IjiSj. However, this single- molecule 
imaging approach is restricted to the study of polymers 
which are larger in size than the diffraction limit of the 
imaging optics. Fluorescence correlation spectroscopy 
(FCS) is a technique that can purportedly resolve dy- 
namics below the diffraction limit. In a recent paper, 
Shusterman et al. [6j introduced an analysis method 
for FCS data, which they used to analyze the nanoscale 
conformational dynamics of single- and double-stranded 
DNA, reporting Zimm instead of Rouse dynamics for the 
studied molecules. This finding has spurred considerable 
theoretical work in an effort to explain this unexpected 
result [IHQ]. Meanwhile, the analysis method of Shus- 
terman et al. has found widespread application in FCS- 
based studies of polymer dynamics Ref. [ToHll] . 

Here, we present a fundamental analysis of the capabil- 
ity of FCS to elucidate polymer dynamics on the nanome- 
ter scale. If FCS is indeed capable of measuring the con- 
formational dynamics of small molecules close to and be- 
yond the diffraction limit of the used optics, this would 
be of tremendous interest. Not only could one resolve 
the question of whether Rouse or Zimm dynamics is the 
appropriate model for the dynamics of disordered poly- 
mers, but one could also experimentally study problems 
such as protein folding and unfolding, or the nanoscale 
diffusion of molecules within lipid membranes. 

Thus, the core question is how well FCS can elucidate 
nanoscale conformational/diffusional dynamics. Let us 
start with a simple one-dimensional system, namely the 
one-dimensional free diffusion, with diffusion coefficient 
Dq^ of a molecule along the x-axis within the confined 
region —a<x<a, where the confinement itself diffuses 
with a much slower diffusion coefficient D through the 
focus of a FCS system. For the diffusion within the con- 
fined region. Green's function is given in a standard way 



as an expansion over eigenfunctions 
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where the eigenfunctions (/>j(x) are non-zero only for 
\x\ < a, and where they are ^2j(^) = cos{J7rx/a) for even 
indices, and (j)2j-\-i{x) = sin((j + l/2)7rx/a) for odd in- 
dices. The characteristic frequencies are Uj = (J7r/2a)^. 
Let us further assume that fluorescence is excited by fo- 
cusing a plane wave through a lens with numerical aper- 
ture NA into a medium with refractive index n. The 
resulting intensity distribution is given by 
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where /Cmax = 27rNA/A. For 500 nm wavelength light, a 
numerical aperture of 1.2, and a refractive index of 1.33 
(water), one can estimate the focus width and thus spa- 
tial resolution of a corresponding scanning microscope to 
be of ca. 200 nm. The best fit of a Gaussian distribu- 
tion to the actual intensity distribution, Eq. Q, is found 
to have a Gaussian width of 2<j = 152 nm. Assuming 
that the excited fluorescence is detected uniformly over 
the whole excitation region, then U{x) is directly pro- 
portional to the molecule detection function (MDF) of 
the system, which gives the probability density of de- 
tecting a photon from a molecule at position x. The 
autocorrelation function (ACF) of an FCS measurement 
is then given as the multiple integral over the product 
of the probability to detect a photon from a molecule at 
some initial position, the probability density that it dif- 
fuses from this position to a final position within time t 
(given by Green's function), and the probability to detect 
a photon from a molecule at this final position. Thus, the 
short-time behavior of the ACF can be written as 
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where y is the center position of the confined region with 
respect to the excitation focus. The integrations over x 
and x^ run from —a to a, and the integration over y from 
— oo to oo. In Eq. ([s]), we have neglected ah constant 
pre-factors related to concentration, overall detection ef- 
ficiency, etc. "Short-time behavior" means that we con- 
sider here only a time range where the fast confined diffu- 
sion dominates the ACF, before the slow diffusion of the 
confined region itself plays any role. Integration over y 
takes into account that the confined region can be at any 
position with respect to the excitation focus. Eq. ([3| can 
be rewritten by switching to Fourier space, representing 
all functions by their Fourier transforms, e.g. 
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After carrying out several integrations, one finds the com- 
pact expression 



g{t) = ^aj exp{-ujjDot) 
where the amplitudes aj are given by 



(5) 
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A: (1/100 nm) 

Fourier power spectra of the MDF and the eigenfunc- 
'j entering the Green's function for confined diffusion. 



0, 1, ..., represent the power spectra |0j(A;)| 



for a = 50 nm. The triangular shaded region is the power 
spectrum |[/(A:)|^, and the yellow shaded region that of the 
Gaussian approximation of U{x). Vertical lines delimit the 
finite support of the optical transfer function. The position 
of the maxima of the functions \(j)j{k)\^ scales with the in- 
verse value of a. Thus, the smaller the confinement region, 
the more the maxima will be shifted away from the Fourier 
modes that can be probed by the MDF. 



The core idea of Shusterman et al. [6] is to compare 
a measured ACF with the model ACF for free-diffusion 
in the Gaussian approximation of the MDF. For the one- 
dimensional case considered here, such a model ACF is 
given, up to a constant pre-factor, by the simple expres- 
sion 
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Eqs. ([5| and (|6| are particularly suited for analyzing the 
impact of the optical resolution of the measurement setup 
on the temporal behavior of an ACF. As can be seen, the 
amplitude for the jth exponential function in Eq. ([5| is 
given by the weighted integral over the power spectrum 
of the jth eigenfunction (pj , where the weight function is 
the power spectrum of the MDF. However, for increas- 
ing numbers of j, the power spectrum of (j)j is shifted to 
higher and higher spatial frequencies /c, whereas the filter 
U{k) is non-zero only for k < /Cmax- A visualization of 
this fact is presented in Fig. [l} showing the power spectra 
for the first few eigenfunctions together with the power 
spectrum of the MDF (both for the exact MDF as well for 
its Gaussian approximation) . The calculations were done 
for a = 50 nm, thus confining diffusion to a region of 100 
nm. As can be seen, there is substantial overlap between 
\U{k)\'^ and \(j)j{k)\'^ for the first three eigenfunctions only 
- all the information contained in the higher order eigen- 
functions which correspond to faster time scales (larger 
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will not be refiected in the ACF. 
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where {x{t)'^) = 2Dot is the mean square displacement 
of a molecule within time t, and a the variance of the 
Gaussian that is used for approximating the MDF. Under 
the assumption that this form of the ACF will also be 
valid for more complex diffusion behavior, Eq. ^ can be 
inverted to yield the alleged relationship 



{x{tf 



gHt) 



(8) 



For free diffusion, the right hand side of Eq. ([8| is a linear 
function of time. For confined diffusion, one would still 
hope to see the linear behavior on short time scales. Let 
us define a quantity a{t) as 
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which defines the local exponent with which (x(t)^) in- 
creases with time t. This function is time-independent 
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FIG. 2: Behavior of the function a{t) for different values of 
the confinement parameter a and for Dq — 1 /xm^/time unit. 
The inset shows the different ACFs for the selected values of 



and equal to one for free diffusion, and one also expects 
that it will be close to one on sufficiently short time scales 
for confined diffusion. Fig. [2] shows the behavior of a{t) 
for several selected values of the confinement parame- 
ter a. Indeed, a{t) converges to unity in the limit of 
zero time, but this is true for any well-behaved function 
g(t)\ However, not much can be learned from a(t) for 
intermediate time values. This is due to the fact that 
with increasing confinement (decreasing value of a) the 
correlation function captures less and less of the diffu- 
sion dynamics due to the finite support of the MDF in 
Fourier space. It is important to notice that we consider 
here only the limiting case of negligible diffusion of the 
confined area as compared to the diffusion of a molecule 
within the confined area. At large times, the diffusion 
of the whole confinement region will start to show up in 
the ACF, thus leading again to values of a{t) close to 
unity. In the intermediate region, however, the curves 
in Fig. [2] show that one can extract any fractional power 
dependence of {x{t)'^) on time t. 

To consider a more realistic scenario, and to prove that 
the behavior of a{t) is not an artifact of the simple one- 
dimensional model system considered so far, let us study 
the rapid diffusion of a particle in three dimensions within 
an isotropic harmonic potential, VCR) = ksTR^ /2i<i^ 
where ks is Boltzmann's constant, T is temperature, and 
/t: is a constant having units of length squared which de- 
termines the steepness of the potential. The coordinate 
vector R refers to the center position Yc of the potential, 
R = r — Tc, whereas the center position Yc is assumed 
to diffuse freely through space (slow diffusion). Denot- 
ing the diffusion coefficient for the rapid diffusion again 
by Do, Green's function for the diffusion within the har- 
monic potential is given by 
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where we have used the abbreviation s = e~^o^/'^. Fur- 
thermore, for the slow diffusion of the center position Yc 
with diffusion constant D < Dq^ we have the Green's 
function for the free diffusion equation 
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In most FCS publications, and also in Ref.[6l [TT], the 
three-dimensional MDF, t/(r), of a confocal fluorescence 
microscope is approximated by an axi-symmetric Gaus- 
sian distribution with half axes a and (a. For calculating 
the final ACF, we still need the equilibrium probability 
Po(Ro) of finding a molecule at some initial position Rq 
within the harmonic potential, which is found by letting 
time t in Eq. ( p!Q| ) approach infinity. The ACF is then 
equal to the multiple integral of the product of the initial 
probability ^0(^0)7 the probability density G(R, Ro,t) 
that the molecule moves from initial position Rq to posi- 
tion R within the harmonic potential during time t, the 
probability density Gc(rc — r^o,^) that the center posi- 
tion diffuses from rc,o to position Yc within the same time 
t, and the probabilities U{yo) and U{y) to detect a pho- 
ton at the initial and final position. After carrying out all 
twelve integrations over all possible initial and final posi- 
tions of molecule and potential center, one finally arrives 
at the analytic expression (up to a constant pre-factor) 
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As an example, the correlation function is calculated for 
the following parameters: the shape of the MDF is de- 
termined by cr = 150 nm and C = 5, the fast diffusion 
coefficient within the potential is I^o = 1 /im^/time unit, 
the diffusion coefficient D of the diffusion of the cen- 
ter position of the harmonic potential is set equal to 
0.5 /im^/time unit, and hz is chosen to be (150 nm)^. Sim- 
ilar to the one-dimensional case, one implicitly defines an 
effective mean-square displacement function (x^) by the 
relation ^ 
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Fig. [3] shows the behavior of (x^) as a function of time, 
together with the corresponding shape of the ACF 
(inset). The broken line represents a linear fit of \nt 
to the curve of In(x^) at the point of minimum slope. 
For the chosen parameters, this slope value is ca. 2/3, 
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Ref. [6] can produce any power-law behavior with Ut- 
tle connection to the actual underlying physical process. 
However, as was shown when analyzing the interplay be- 
tween the Fourier spectra of the eigenfunctions entering 
the Green's function of the diffusion equation and the 
Fourier spectrum of the MDF (Fig. fl]), an FCS exper- 
iment cannot capture much of the dynamics on length 
scales that are outside of the finite Fourier support of the 
MDF of the optical system. Thus, our analysis highlights 
the general limitation of FCS for elucidating processes on 
length scales below the resolution limit of the optics used 
in an FCS experiment. 



FIG. 3: Mean-square displacement (x^) as a function of time 
determined by inverting Eq. (14) for the parameter values 
given in the text. Inset shows the corresponding ACF as 
calculated from eqs |12| and |13| The dashed line is a linear fit 
of Int to In(x^) at the point of minimum slope. 
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FIG. 4: Minimum value of a{f) as a function of the ratio 
D/Dq. Shown are also the expected values of a for Zimm 
and Rouse dynamics. 



mimicking a diffusion process originating from a Zimm 
dynamics scenario. Let us define again a characteristic 
function a{t) as in Eq. ([9|. Generally, a{t) will approach 
unity for very small and very large values of time t, but 
will have a minimum at intermediate time values. Fig. |4] 
shows the minimum value of a(t) as a function of the 
ratio D/Dq. As can be seen, depending on this value, 
the function {x^) can show intermediate power-law 
behavior with any exponent, although the underlying 
process is simple diffusion within a harmonic potential. 

In summary, FCS measurements are governed by the 
same spatial resolution limit as any linear optical imag- 
ing system, and the smaller the length scale of a diffu- 
sion process the less information FCS can capture about 
it. Moreover, defining and analyzing an effective mean- 
square displacement function {x^) as first proposed in 
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